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ABSTRACT 

The use of Computer tools to solve calculations in any area of study is 
essential to reduce the time spent performing them, leading to the investment of time 
in the analysis of the problem for a number "n" of Solutions to the proposed problem, 
being able to perform iterations múltiple until reaching an optimal result. The 
objective of this article is to use one of the simplest Computer tools, such as the 
Microsoft Excel program, to create a code in Visual Basic that allows the resolution 
of fíat structures, subject to shear forces, axial forces and bending moments, thus 
calculating the reactions and displacements in the nodes of the structure. The macro 
was based on the analysis of structures with the directstiffness method, where codes 
were created to obtain the stiffness matrix of each element, the global stiffness matrix 
of the structure assembled from the stiffness matrix of each element, the matrix of 
embedding forces, the matrix of nodal forces, and the matrix of displacements of 
each node. Each of the matrices are displayed in a spreadsheet, so that for academic 
purposes, the user can view each of the matrices. 

Keywords: Excel Macros, visual basic, planar structures, stiffness matrix. 

USO DE MACROS EN MICROSOFT EXCEL PARA 
ANALIZAR ESTRUCTURAS PLANAS. 

RESUMEN 

El uso de herramientas informáticas para resolver cálculos en cualquier área 
de estudio es esencial para reducir el tiempo dedicado a realizarlas, dando lugar a 
la inversión de tiempo en el análisis del problema para un número "n" de soluciones 
al problema propuesto, pudiendo realizar iteraciones múltiples hasta alcanzar un 
resultado óptimo. El objetivo de este artículo es utilizar una de las herramientas 
informáticas más simples, como el programa Microsoft Excel, para crear un código 
en Visual Basic que permita la resolución de estructuras planas, sujetas a fuerzas 
de corte, fuerzas axiales y momentos flectores, calculando así las reacciones y 
desplazamientos en los nodos de la estructura. La macro se basó en el análisis de 
estructuras con el método de rigidez directa, donde se crearon códigos para obtener 
la matriz de rigidez de cada elemento, la matriz de rigidez global de la estructura 
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ensamblada a partir de la matriz de rigidez de cada elemento, la matriz de fuerzas 
de empotramiento, la matriz de fuerzas nodales, y la matriz de desplazamientos de 
cada nodo. Cada una de las matrices se muestran en una hoja de cálculo, de modo 
que para fines académicos, el usuario podrá visualizar cada una de las matrices. 

Palabras clave: Macros de Excel, Visual Basic, estructuras planas, matriz de 
rigidez. 

1. INTRODUCCION. 

La resolución de estructuras planas por métodos matriciales, resulta ser un 
trabajo de acuerdo con el número de elementos (barras y nodos) presentes en la 
estructura, lo que representa un gran número de cálculos por elemento que, al 
trabajar manualmente, podrían incurrir en errores, como como el uso de pocos 
decimales, transponer valores, transponer índices matrices que alterarían el orden 
de los coeficientes de la matriz de la estructura y por lo tanto el cálculo erróneo de 
los valores de rigidez. Por lo tanto, este trabajo de investigación se basa en la 
programación de macros en Microsoft Excel para el cálculo de estructuras planas, 
con un enfoque estrictamente académico basado en métodos matriciales, con el fin 
de que el estudiante dedique más tiempo a analizar el problema que a resolver el 
problema. Sin quitar la importancia de conocer y realizar los cálculos manualmente. 

Se sabe que los programas existentes utilizan métodos de cálculo matricial 
para resolver estructuras. Sin embargo, las matrices generadas por cada elemento 
no pueden ser visualizadas por el usuario, solo los resultados obtenidos de la 
resolución final, por lo que académicamente es difícil para el estudiante tener un 
punto de comparación para verificar sus ejercicios, y esto se agrega a eso Los 
programas que están en el mercado tienen licencias de uso caras y temporales, que 
probablemente no todos los estudiantes poseen, y además se requiere que el 
estudiante tenga conocimiento del uso de los programas para poder manipularlo. 

2. BASES TEÓRICAS. 

La macro se desarrolló basándose en los métodos de análisis matricial de 
estructuras, específicamente con el método de rigidez directa. La formulación 
general de este método se muestra en la ecuación 1, de acuerdo con (Uribe 
Escamilla, 2000), (Weaver & Gere, 1990), (Aguiar, 2014) y (Mueller, 2016) que se 
utiliza para calcular las fuerzas de reacción en los apoyos. 

{F}=[K]*{D}+m (1) 

Donde: 

{F} es el vector de fuerzas nodales (reacciones y fuerzas puntuales). 

{K} es la matriz de rigidez. 

{D} es el vector de desplazamientos nodales. 

{Qe} es el vector de fuerzas de empotramiento. 

El cálculo de los desplazamientos nodales se realiza de acuerdo a la ecuación 2. 
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{£>} = [ K - 1 ] * {F - Qf) (2) 

Donde [ K -1 ] es la matriz de flexibilidad de la estructura. 

El método de cálculo en el que se basa la programación de la macro es el método 
de los desplazamientos. 

2.1. Solución general por el método de los desplazamientos. 

Partiendo del enfoque anterior, se realiza una matriz basada en fuerzas y 
desplazamientos conocidos y desconocidos, guardando la relación de que con las 
fuerzas conocidas, los desplazamientos son desconocidos y viceversa. El orden se 
muestra en la ecuación 3 según (Uribe Escamilla, 2000) y (Weaver & Gere, 1990). 

í0n\(K nn Kna\jDn\ (3) 

\Qa) v Kan Kaa' \Da) 

Donde: 

Q n es el vector de fuerzas nodales conocidas. 

Qa es el vector de fuerzas nodales desconocidas (reacciones). 

Knn es la submatriz de rigidez utilizada para el cálculo de los desplazamientos 
nodales. 

Kna es la submatriz de rigidez utilizada para el cálculo de las reacciones. ( Kna = 
Kan). 

Kaa es la submatriz de rigidez utilizada para el cálculo de las reacciones debidas a 
desplazamientos en los apoyos. 

Dn es el vector de desplazamientos nodales conocidos. 

Da es el vector de desplazamientos nodales desconocidos. 

3. MATRIZ DE RIGIDEZ PARA UN ELEMENTO ORIENTADO 
ARBITRARIAMENTE. 

En general, un elemento estructural 2D está sujeto a fuerzas axiales, 
fuerzas de corte y momentos de flexión. Los ejes originales del elemento son ejes 
locales, mientras que los ejes de toda la estructura son ejes globales, por lo que es 
necesario transformar los ejes locales en ejes globales, como se muestra en la 
Figura 1. 
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Figura 1 . Caso general de un elemento estructural 2D orientado 
arbitrariamente. (López, 2011) 

La matriz de transformación (T) entre ejes locales y globales se muestra en 
la ecuación 4. 

eos <fvc sen(¡K 0 0 0 0 

- sent/K eos <jx 0 0 0 0 

0 0 1 0 0 0 

(41 

0 0 0 eos (¡k sen<¡K 0 v ' 

0 0 0 - sentjfx costfvc 0 

0 0 0 0 0 1 


La matriz de rigidez de un elemento estructural orientado arbitrariamente, 
puede obtenerse de una barra horizontal, sometida a fuerzas de corte y momentos 
flectores, cuando se generan los siguientes desplazamientos y rotaciones en sus 
extremos. 


3.1 Desplazamiento vertical unitario en el nodo i. 


Las reacciones en los nodos de la viga estudiada se pueden obtener 
aplicando las relaciones de Mohr (Beer, Johnston, DeWolf, & Mazurek, 2010), 
(Hibbeler, 2012), (Popov, 2000) como se muestra en las ecuaciones 5, 6, 7 y 8. 


1 f x 

V(x) = — J ~ wdx 

M(x) = Y¡j v(x)dx 
1 r° x 

0(x) = — M(x)dx 
J o 
l r x 

y(x) = —J 0(x)dx 


(5) 

( 6 ) 

(7) 

( 8 ) 


Resolviendo las ecuaciones anteriores, se obtienen las ecuaciones 9, 10, 
11 y 12 (conw = 0). 
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V(x) = 


Cí 

T 

i 


M(x) = — (Clx + C2) 

1 

0(x) = — (Clx 2 + C2x + C 3) 

El 

1 

y (x) = — (Clx 3 + C2x 2 + C3x + C4) 
ti 


(9) 

( 10 ) 

( 11 ) 

( 12 ) 


Las ecuaciones 9,10,11 y 12, se resuelven utilizando las siguientes condiciones de 
frontera: 

Tabla 1. Condiciones de frontera para la resolución de las ecuaciones. 


X 

y 

0 

0 

vi = 1 

0 

L 

0 

0 


Las constantes C1, C2, C3 y C4 poseen los siguientes valores: 


El 


El 


C 1 = 12 * — ;C2 = 3 = 0; C4 = El 

Sustituyendo los valores de C1 y C2 en las ecuaciones 9 y 10, se obtienen las 
ecuaciones 13 y 14. 


, . El 

F(x) = 12 * ^ 

. . El El 

M(x) = 12 * -jjx — 6 * -jj 


(13) 

(14) 


Resolviendo las ecuaciones 13 y 14, para x = 0 and x =L, se obtienen las fuerzas 
en los nodos, como se muestra en la Figura 2. 



ffC 


Figura 2. Elemento horizontal sometido a fuerzas de corte y momentos flectores 
debido a desplazamiento vertical unitario en el nodo i. (López, 2011) 

3.2 Desplazamiento vertical unitario en el nodo j. 

Aplicando el mismo procedimiento anterior, y variando las condiciones de 
contorno para un desplazamiento vertical unitario en el nodo j, se obtienen los 
valores de las reacciones en los nodos como se muestra en la Figura 3. 
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Figura 3. Elemento horizontal sometido a fuerzas de corte y momentos flectores 
debido a desplazamiento vertical unitario en el nodo j. (López, 2011). 

3.3 Rotación unitaria en el nodo i. 

Aplicando el mismo procedimiento anterior, y variando las condiciones de 
contorno para una rotación unitaria en el nodo j, se obtienen los valores de las 
reacciones en los nodos como se muestra en la Figura 4. 


60 

L 2 


t 



Figura 4. Elemento horizontal sometido a fuerzas de corte y momentos flectores 
debido a rotación unitaria en el nodo i. 

3.4 Rotación unitaria en el nodo j. 

Aplicando el mismo procedimiento anterior, y variando las condiciones de 
contorno para una rotación unitaria en el nodo j, se obtienen los valores de las 
reacciones en los nodos como se muestra en la Figura 5. 


6EI 

L 2 


t 



Figura 5. Elemento horizontal sometido a fuerzas de corte y momentos 
flectores debido a rotación unitaria en el nodo j. 
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3.5 Desplazamiento horizontal unitario en el nodo i y en el nodo j. 

Basándose en la Ley de Hooke, la tensión axial para un elemento estructural 
sometido a fuerzas de compresión o tracción, se muestra en la ecuación 15. 


EA 

P = T 


(15) 


La matriz de rigidez para un elemento estructural, en ejes locales, sometida a 
fuerzas de corte, fuerzas axiales y momentos de flexión, se muestra en la ecuación 
16, según (Uribe Escamilla, 2000), (Weaver & Gere, 1990), (Gavin, 2014), y (Kwon 
& Bang, 1996). 


AE 

L 

o 

o 

_AE 

L 

o 

o 


o 

12 El 
6EI 

H 

0 


0- 


L 

0 


L 


6 El 

0 

12 El 

6 El 

~zT 

~¡r 

n 

4EI 

0 

6 El 

2 El 


AE 

L 


L 

0 


L 

0 


12 El 

6 El 

0 

12 El 

6 El 

ú 


ú 

L 2 

6 El 

2 El 

0 

6EI 

4 El 

IT 

L 

n 

L 


(16) 


Para obtener la matriz de rigidez de un elemento orientado arbitrariamente, 
es necesario convertir los ejes locales en ejes globales, aplicando la ecuación 17, 
según (Uribe Escamilla, 2000), (Weaver & Gere, 1990), (Gavin, 2014), and (Kwon & 
Bang, 1996). 


K = T T * K * T 


(17) 


Una vez aplicada la ecuación 17, la matriz de rigidez para un elemento 
estructural orientado arbitrariamente se muestra en la ecuación 18. 


Tae 

I2E!]. 

6EI 

[AE., 12EI ,1 

f AE 

12EI1 

ir 

~írr 

“¡r 


Ít 

~~¡r\ 


[AE 

jml 

Lr 

L 3 J 


6FA 


AE 2 12EI ,i 
—ir +— j-tr 
L L 3 


GEI, 






{ 


AE ]_2ET 
L " Ü 

6E1 




í AE > 12EI , 
[ L H L 3 


ÓEL 



[AE IZEll 

- 

IT' ~\ 

4 El 

6ET 

L 

L 2 ^ 

6EÍ AE o 12 El 

^ ___X + —s-j 

L L L 3 

1 6E 1 I 

'AE E2EI1 

J ~¡r k 1 

. L " L 3 J 


2EI 

L 


6E1 


n 12E1 ,1 6EL 


JL 

L 2 A 


12El i [AE 12EII 


AE 2 12E1* j 

— u + —*—A 

L ^ L 


6EL 

~TT X 


2F1 

l 


6 El 


6E1. 

“17 a 


4EI 

L 


( 18 ) 
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Donde: 

A es el área transversal del elemento en estudio. 

E es el módulo de elasticidad del material del elemento estructural.. 

/ es el momento de inercia de la sección transversal. 

L es la longitud del elemento estructural. 

Y es el cos0x. 

H es el sin 0x. 

4 DESARROLLO DE LA MACRO. 

La programación de la hoja de cálculo se realizó con macros en Microsoft 
Excel. Una macro en Microsoft Excel se define como la interacción entre Microsoft 
Excel y Visual Basic en la misma hoja de cálculo, pudiendo realizar operaciones 
complejas en la hoja de cálculo de forma rápida y sencilla. 

El programa de ensamblaje consta de 3 hojas de Excel, donde la primera 
se denomina "Matriz de rigidez", y en ella se introducen los datos de geometría de 
la estructura en estudio para generar la matriz de rigidez de cada elemento y luego 
ensamblar la matriz de rigidez del conjunto completo, estructura; otra hoja llamada 
"Matriz de Fuerzas", donde se calculan las fuerzas de incrustación de los elementos, 
se introducen las fuerzas nodales y los desplazamientos nodales conocidos; una 
hoja llamada "Cálculos”, donde se realizan las operaciones matriciales necesarias 
para calcular las reacciones y los desplazamientos. 

4.1 Programación de la hoja "matriz de rigidez". 

En la hoja de entrada de datos, se solicitarán los datos de geometría y 
topología de la estructura en estudio, y los datos de entrada serán los indicados en 
Tabla 2. 


En esta etapa, la macro solo está programada para leer los datos de 
entrada, la macro es responsable de leer estos datos con las coordenadas de las 
filas en números y las columnas en números, por ejemplo, la celda Al en Excel, es 
la celda (1,1) en el código de Visual Basic, porque las columnas no se leen como 
letras sino en números. El comando en Visual Basic para leer y asignar valores de 
celdas es "Celdas (fila, columna)" (Shepherd, 2004), como se muestra en el 
apéndice 1, donde "i" es el número de elementos de la estructura, que se lee 
automáticamente del número de filas de los datos de entrada. Los valores de la 
matriz de rigidez para cada elemento tienen un código en Visual Basic como se 
muestra en el apéndice 1, donde x y w, son valores de filas y columnas para impri mir 
los valores calculados en Excel. Para el ensamblaje de la matriz de rigidez global, 
se realizó un ciclo doble "For-Next" (Mueller, 2016), para buscar y sumar 
algebraicamente los valores de rigidez para nodos repetidos, como se muestra en 
el apéndice 1. 
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Para calcular la matriz de rigidez de cada elemento, se debe presionar un 
botón denominado "Crear matriz de rigidez de cada elemento", para activar la 
macro, como se muestra en la Figura 6. 


Sección 

Base (m) 

Altura (m) 

¿A/iea (m 2 ) 

Inercia en x 
(m 4 ) 

0.3 

0,3 

0,09 

0.000675 

0.3 

0.4 

0,12 

0.001-6 














NQ de nodos 

N? de b arras 

I v 1 

2 1 


Crear M. 


(riz de Rigidez 


de cada elemento 


Figura 6. Botón para crear la matriz de rigidez de cada elemento. Fuente: Autores. 


Una vez que se hayan creado las matrices de rigidez para cada elemento, 
la macro de ensamblaje debe activarse utilizando otro botón, denominado "Crear 
matriz de rigidez global", como se muestra en la Figura 7. 


theta2 

O 

31476093,75 

83936250 

O 

-31476093,75 

167872500 



Crear Matriz 
Global 




u2 

Q2x 

223830000,00 

M2 

0,00 

Qlx 

-223830000,00 

Qiy 

0,00 

MI 

0,00 

Q2y 

0,00 


Figura 7. Botón para ensamblar la matriz de rigidez global. Fuente: Autores. 

Para probar la precisión de la macro, un ejemplo del libro del autor (Uribe 
Escamilla, 2000), ejemplo 11.22 (ver Figura 8), en la página 520, será resuelto. 


120 kN 



-pf- 


4 m 


-y 


Figura 8. Ejemplo 11.22, página 520. (Uribe Escamilla, 2000). 

La numeración de los nodos para este ejercicio, se muestra en la Figura 

9 ¡ Error! No se encuentra el origen de la referencia.. 
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Figura 9. Numeración de nodos para datos de entrada (Uribe Escamilla, 

2000 ). 

Los datos requeridos para el análisis, se muestran en la Tabla 2.¡Error! No 
se encuentra el origen de la referencia. 

Tabla 2. Datos de entrada de la estructura. Fuente: Autores 


Elemento 

Nodo 

i 

Nodo 

i 

Longitud 

(m) 

6(°) = 

Módulo de 
elasticidad 
(Kg/m 2 ) 

Base 

(m) 

Altura 

(m) 

Área 

(m 2 ) 

Inercia en 
x (m 4 ) 

1 

1 

2 

3 

0 

19000000 

0,3 

0,3 

0,09 

0,000675 

2 

3 

1 

5 

53,13 

19000000 

0,3 

0,4 

0,12 

0,0016 


Las matrices de rigidez de cada elemento se muestran en la Tabla 3. 
Tabla 3. Matriz de rigidez de cada elemento (kN/m).Fuente: Autores. 


Elemento 1 

ul 

vi 

thetal 

u2 

v2 

theta2 

Qlx 

570000 

0 

0 

-570000 

0 

0 

Qly 

0 

5700 

8550 

0 

-5700 

8550 

MI 

0 

8550 

17100 

0 

-8550 

8550 

Q2x 

-570000 

0 

0 

570000 

0 

0 

Q2y 

0 

-5700 

-8550 

0 

5700 

-8550 

M2 

0 

8550 

8550 

0 

-8550 

17100 

Elemento 2 

u3 

v3 

theta3 

ul 

vi 

thetal 

Q3x 

166027,78 

217479,17 

-5836,8 

-166027,8 

-217479,17 

-5836,8 

Q3y 

217479,17 

292890,62 

4377,6 

-217479,2 

-292890,62 

4377,6 

M3 

-5836,8 

4377,6 

24320 

5836,8 

-4377,6 

12160 

Qlx 

-166027,8 

-217479,2 

5836,8 

166027,78 

217479,17 

5836,8 

Qly 

-217479,2 

-292890,6 

-4377,6 

217479,17 

292890,62 

-4377,6 

MI 

-5836,8 

4377,6 

12160 

5836,8 

-4377,6 

24320 


La matriz de rigidez de cada elemento según (Uribe Escamilla, 2000), se 
muestran en la Figura 10. 
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U| 


e» 

v t 

«i 

V ! 

e, 

U] 

v i 

0, 


570000 

0 

0 

lili 


Ififiüio 

217480 

- 5840 

-166030 

-217480 

-5840 

y 3 

0 

5700 

8550 

1 l l 


217480 

292890 

4380 

-217480 

-292890 

4380 


0 

8550 

17100 

l l í 


-5840 

4380 

24320 

5840 

-4380 

12160 

% 

-570000 

0 

0 

1 i 1 


-166030 

”217480 

5840 

166030 

217480 

5840 

u, 

0 

-5700 

-H550 

1 1 1 


“217480 

-292890 

-4380 

217480 

292890 

-4380 

v i 

O 

8550 

8550 

1 1 1J 


-5840 

4380 

12160 

5840 

-4380 

24320 



Figura 10. Matriz de rigidez de cada elemento (kN/m). (Uribe Escamilla, 

2000 ) 

Obsérvese que las matrices obtenidas por el macro (Tabla 3) son iguales a 
las calculadas por el libro fuente (Uribe Escamilla, 2000) (Figura 10). 

La matriz global debe ser ensamblada, y manualmente, el proceso podría 
ser largo debido a la suma de los factores de rigidez para nodos repetidos. La macro 
tiene un "ciclo de buscador y acumulador", que funciona buscando los índices en 
cada matriz de rigidez de los elementos y sumando los valores de los índices 
iguales. Al presionar el botón "Crear matriz de rigidez global", la macro calcula la 
matriz global como se muestra en la Tabla 4¡ Error! No se encuentra el origen de 
la referencia.. Los desplazamientos en el nodo 1 son desconocidos, y de acuerdo 
a (Uribe Escamilla, 2000), (Weaver & Gere, 1990), (Mueller, 2016), (Dean), estos 
índices deben organizarse en primer lugar, y luego se organizan los valores de los 
desplazamientos conocidos que corresponden a los de los apoyos en los nodos 2 y 
3. 


Tabla 4. Matriz global de la estructura ensamblada. (López, 2011). 



ul 

Vi 

thetal 

u2 

v2 

theta2 

u3 

v3 

theta3 

Qlx 

736027,78 

217479,17 

5836,80 

-570000,00 

0,00 

0,00 

- 166027,78 

- 217479,17 

5836,80 

Qiy 

217479,17 

298590,62 

4172,40 

0,00 

- 5700,00 

8550,00 

- 217479,17 

- 292890,62 

- 4377,60 

MI 

5836,80 

4172,40 

41420,00 

0,00 

- 8550,00 

8550,00 

- 5836,80 

4377,60 

12160,00 

Q2x 

- 570000,00 

0,00 

0,00 

570000,00 

0,00 

0,00 

0,00 

0,00 

0,00 

Q2y 

0,00 

- 5700,00 

- 8550,00 

0,00 

5700,00 

- 8550,00 

0,00 

0,00 

0,00 

M2 

0,00 

8550,00 

8550,00 

0,00 

- 8550,00 

17100,00 

0,00 

0,00 

0,00 

Q3x 

- 166027,78 

- 217479,17 

- 5836,80 

0,00 

0,00 

0,00 

166027,78 

217479,17 

- 5836,80 

Q3y 

- 217479,17 

- 292890,62 

4377,60 

0,00 

0,00 

0,00 

217479,17 

292890,62 

4377,60 

M3 

5836,80 

- 4377,60 

12160,00 

0,00 

0,00 

0,00 

- 5836,80 

4377,60 

24320,00 


La submatriz que contiene los subíndices referidos al nodo 1, es la matriz 
Knn, como se muestra en la Figura 11, la cual coincide con la matriz creada por la 
macro (las primeras 3 columnas y 3 primeras filas de la Tabla 4¡ Error! No se 
encuentra el origen de la referencia.). 

736030 217480 5840' 

217480 298590 4170 

5840 4170 41420 

Figura 11. Matriz Knn. (Uribe Escamilla, 2000). 


La matriz Kna se muestra en la Figura 12. 
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-570000 

0 

0 ’ 


'-166030 - 217480 

- 5840 ' 

0 

-5700 -8550 


-217480 -292890 

4380 

0 

8550 

8550 


5840 -4380 

12160 


Figura 12. Matriz Kna. (Uribe Escamilla, 2000). 

Se observa que todas las submatrices calculadas con la macro son muy 
cercanas a las matrices generadas en el ejercicio del libro de (Uribe Escamilla, 
2000 ). 


4.2 Programación de la hoja “Matriz de fuerzas”. 

En esta hoja, el usuario debe ingresar las fuerzas aplicadas en los nodos y 
elementos. La hoja reconoce automáticamente el número de elementos y nodos, y 
solo pregunta por el tipo de fuerzas aplicadas, como se muestra en la Figura 13 y 
Figura 14. 


1 D 1 

KB 


| f] | 

Iff T i 


t . 


í^nlllw, 

1 ■■■■■■■■■ 4 

A « 4 • 


_i 

4- « ♦ * l , 

— ♦- fe -l-l 


jU-^-4- > -L 


Figura 13. Tipos de fuerzas en elementos estructurales. Fuente: Autores. 




. r 

4 1 4 i- - 

M m -i- i .| 

~ 1 f ° - b 


Figura 14. Tipos de fuerzas en elementos estructurales. Fuente: Autores. 


Una vez elegido el tipo de fuerzas aplicadas a los elementos, el valor de la 
fuerza debe ingresarse en las celdas correspondientes, como se muestra desde la 
Tabla 5 y Tabla 6. 

Tabla 5. Datos de entrada para la generación de fuerzas de 
empotramiento con la macro. Fuente: Autores. 


Cuadrante 

Elemento 

Nodo 

i 

Nodo 

j 

L (m) 

e(°) = 

Solicitación 

(Kg/m) 

c (m) 

L/2 (m) 

a (m) 

b (m) 

4 

1 

1 

2 

3 

0 

6 

120 



1,5 

1,5 

5 

2 

3 

1 

5 

53,13 

5 

20 


2,5 
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Tabla 6. Fuerzas de empotramiento calculadas con la macro. Fuente: 

Autores. 


Reacciones de empotramiento 
(coordenadas locales) 

Reacciones de empotramiento (coordenadas 
globales) 

Mi 

(Kg.m) 

M j (Kg.m) 

Ri (kg) 

Ry (Kg) 

Rix (kg) 

R/y (Kg) 

Rjx(kg) 

Ftiy(Kg) 

45,00 

-45,00 

60,00 

60,00 

0,00 

60,00 

0,00 

60,00 

41,67 

-41,67 

50,00 

50,00 

0,00 

0,00 

0,00 

0,00 


La fuerzas de empotramiento de acuerdo a (Uribe Escamilla, 2000), se 
muestran en la Figura 15 and Figura 16. 



Xn = X 21 = 0 

Y¡2 = Y 21 = 120/2 = 60.0 kN t 
M 12 =-M 2 i =120x3/8 = 45.0 kN-in) 


Figura 15. Fuerzas de empotramiento en el elemento 1. Fuente (Uribe Escamilla, 

2000 ). 



X' = X* = -50.0 X sen <|> = -40.0 kN <- 
= 50.0x cosí). = 30.0 kN t 

M,, =-Mfj =41.67 kN-nO 


Figura 16. Fuerzas de empotramiento en el elemento 2. Fuente (Uribe Escamilla, 

2000 ). 

Una vez calculadas las fuerzas nodales debidas a las fuerzas en los 
elementos a los nodos, es necesario organizar estas fuerzas según el orden 
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(índices) de la matriz global, es decir, las fuerzas en el nodo 1 deben ser las primeras 
y las fuerzas en los nodos 2 y 3, debe ser el último, justo en ese orden. El Código 
se activa al presionar el botón llamado “Generar matriz de fuerzas y 
desplazamientos” (ver Figura 17), y los resultados se muestran en la Tabla 7. 

I 

Generar Matriz de Fuerzas y 
Desplazamientos 


Figura 17. Botón para generar matrices de fuerzas y desplazamientos. Fuente: 

Autores. 

Tabla 7. Matrices de fuerzas de empotramiento calculadas con la macro. Fuente: 

Autores. 


Elemento 1 

Qlx 

0,00 

Qiy 

60,00 

MI 

45,00 

Q2x 

0,00 

Q2y 

60,00 

M2 

-45,00 


Elemento 2 

Q2x 

-40,00 

Q2y 

30,00 

M2 

41,67 

Q3x 

-40,00 

Q3y 

30,00 

M3 

-41,67 


Observe que la matriz que se muestra en la Tabla 7, está organizada para cada 
elemento. La matriz organizada de fuerzas y desplazamientos se calcula mediante 
la macro (ver Tabla 8). Los valores en la matriz de fuerzas indican que el nodo tiene 
fuerzas puntuales en él, mientras que las celdas que no son valores indican que las 
fuerzas en el nodo son desconocidas (reacciones). En la matriz de desplazamiento, 
los valores con "ceros" indican que los desplazamientos y las rotaciones en los 
soportes son iguales a cero, mientras que las celdas sin valores indican que los 
desplazamientos son desconocidos. 

Tabla 8. Matriz de fuerzas y desplazamientos organizadas. Fuente: Autores. 


Fuerzas de 
empotramiento 

Fuerzas 

puntuales 

Q2x 

0,00 

100 

Q2y 

-60,00 

0 

M2 

-45,00 

0 

Q3x 

-40,00 


Q3y 

30,00 


M3 

41,67 


Q4x 

0,00 


Q4y 

0,00 


M4 

0,00 



Matriz de 
desplazamientos 

u2 


v2 


theta2 


u3 

0 

v3 

0 

theta3 

0 

u4 

0 

v4 

0 

theta4 

0 
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Con esta convención, la macro detecta que las celdas vacías son incógnitas. 

4.3 Programación de la hoja “Cálculos” 

Básicamente, esta hoja calcula los desplazamientos y reacciones de la 
estructura. Los desplazamientos se calculan utilizando la matriz inversa de rigidez 
(Weaver & Gere, 1990). El inverso de esta matriz, no debe tener el orden completo 
de la matriz de rigidez global, solo tiene el orden del número de desplazamientos 
desconocidos. El código de Visual Basic para calcular la matriz inversa, los 
desplazamientos en los nodos y las reacciones en los soportes, se muestra en el 
Apéndice 1. 

Para activar el código, se debe presionar un botón llamado "Calcular 
reacciones y desplazamientos", como se muestra en la Figura 19. 


Calcuíaze reactions and displacemems 


Figura 18. Botón para calcular y reacciones y desplazamientos. Fuente: Autores. 

Once pushed the button, displacements and reactions are calculated, as shown in 
Tabla 9 and Tabla 10. 


Tabla 9. Matriz inversa. Fuente: Autores. 


Matriz inversa 

1/73179E-06 

-l,25972E-06 

-l,17143E-07 

-l,25972E-06 

4,27011E-06 

-2,52629E-07 

-l,17143E-07 

-2,52629E-07 

2,41849E-05 


Tabla 10. Matriz de resultados. Fuente: Autores. 


Resultados 

0,000254 

ul? 

-0,000371 

vi? 

-0,001085 

thetal? 

-184,80 

Q2x? 

101,39 

Q2y? 

-15,78 

M2? 

4,80 

Q3x? 

78,61 

Q3y? 

-51,75 

M3? 


Los resultados de la estructura según (Uribe Escamilla, 2000), se muestran 

en la Figura 19. 
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u, -0.356 x 10 J m -*• 

v, =—0.560 x 10 ’ m i 
0, = -0.736 x 10 4 rad 


X =X,' 


'-20105 kN *- 


= 

63.83 kN T 

_Mj =Mj,_ 


-5042kN-m 

X =x„' 


' 23.04 kN -+ 

Y, = Y„ 

= 

11617 kN T 

_Mj =Mj,_ 


45.30 kN-rrO 


Figura 19. Resultados de la estructura. (Uribe Escamilla, 2000). 

Obsérvese que la macro convergió a los resultados correctos, coincidiendo 
con la solución plantada por (Uribe Escamilla, 2000). 

5 CONCLUSIONES. 

Como conclusión, se puede establecer que la macro en Excel funciona y se 
ejecuta de manera simple y didáctica para la resolución de estructuras mediante el 
método de rigidez directa y con una programación simple en un programa conocido 
y accesible por cualquier usuario, para que pueda ser usado con fines académicos 
pudiendo observar la generación de todas las matrices inmersas en los cálculos. 

Esta herramienta informática se puede utilizar como complemento de las 
actividades académicas del estudiante y la resolución de los ejercicios, pudiendo 
dedicar más tiempo a analizar el problema que a realizar una gran cantidad de 
cálculos. 
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